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Abstract — The  remote  sensing  group  of  the  University  of 
Hamburg  is  working  in  the  field  of  HF  radar  since  1980.  For 
the  start  three  CODAR  systems  have  been  purchased  from 
NOAA/ERL  (developed  by  D.  Barrick’s  NOAA  group).  Based 
on  16  years  of  experience  a  new  system  called  WEllen  RAdar 
(WERA)  has  been  designed  at  the  University  of  Hamburg  in  1996. 
The  new  design  aims  to  be  as  flexible  as  possible  in  order  to  allow 
easy  adjustment  to  different  requirements,  i.e.  working  frequency, 
spatial  resolution,  and  antenna  configurations.  The  first  part  of 
this  paper  describes  the  technical  solutions  available  to  achieve 
resolution  in  range  and  azimuth.  Modulation  techniques  for  range 
resolution  like  Pulses  and  Frequency  Modulation  (FMCW)  are 
compared,  as  well  as  Direction  Finding  and  Beam  Forming  for 
azimuthal  resolution.  A  short  introduction  to  the  algorithms  is 
given.  The  second  part  discusses  the  hardware  and  software 
components  which  form  a  WERA  and  are  now  commercially 
available.  The  third  part  shows  an  example  of  a  monitoring 
system  bringing  together  HF  radar  remote  sensed  data  and 
numerical  models. 

I.  Introduction 

In  1980,  one  year  after  the  Marine  Remote  Sensing  (MAR- 
SEN)  Experiment  [3],  the  University  of  Hamburg  Remote 
Sensing  Group  started  working  on  HF  radar.  In  close  coope¬ 
ration  to  the  HF  radar  group  at  NOAA1,  at  that  time  lead 
by  D.  Barrick  [1],  three  Coastal  Ocean  Dynamics  Application 
Radars  (CODAR)  systems  have  been  purchased.  The  NOAA- 
CODAR  is  completely  different  from  the  actual  SeaSonde 
which  is  now  delivered  from  D.  Barrick’s  company  “CODAR 
Ocean  Sensors”.  It  uses  a  Continuous  Wave  (CW)  pulsed 
modulation  scheme  for  range  resolution  and  4  receive  antennas 
arranged  in  a  sqare  for  azimuthal  resolution. 

In  1981,  the  first  experiment  with  the  University  of  Hamburg 
CODAR  took  place  on  the  island  of  Sylt,  Germany  [4].  Until 
1983,  several  modifications  have  been  done  to  improve  the 
sensitivity.  The  working  frequency  has  been  increased  from 
25.4  MHz  to  29.85  MHz  in  order  to  reduce  the  impact  from 
radio  interference  due  to  ionospheric  reflections  and  between 
1983  and  1990,  a  shipborne  version  [7]  of  CODAR  has  been 
developed. 
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In  autumn  1983,  the  Canadian  Memorial  University  of 
Newfoundland  in  St.  John’s  organized  an  “International  Work¬ 
shop  on  the  Remote  Sensing  of  Oceanic  Variables  Using  HF 
Groundwave  Radar”.  On  this  workshop  E.  D.  R.  Shearman 
presented  results  of  his  PISCES  HF  radar  system  [13],  which 
used  “Frequency  Modulated  Interrupted  Continuous  Wave” 
(FMICW)  modulation  for  range  resolution  and  a  large  linear 
array  of  receiving  antennas  for  azimuthal  resolution.  Due 
to  the  low  working  frequency  around  9  MHz,  the  working 
range  of  this  system  was  as  large  as  150  km  for  ocean  wave 
measurements  and  even  more  for  currents.  This  presentation 
provided  the  initial  ideas  for  the  University  of  Hamburg 
WEllen  RAdar  (WERA),  however  it  needed  another  12  years 
to  bring  together  the  complete  design  and  the  money  to  build 
a  prototype  WERA  [8].  Meanwhile  a  commercial  version  of 
WERA  is  available. 

II.  Spatial  resolution  of  an  HF  radar 

As  with  every  radar,  spatial  resolution  is  required  to  di¬ 
stinguish  targets  or,  in  case  of  a  radar  for  oceanography, 
different  patches  of  the  sea  surface.  The  following  sections 
describe  modulation  techniques  for  range  resolution  of  an  HF 
radar  as  well  as  antenna  designs  and  algorithms  for  azimuthal 
resolution.  The  most  complete  discussion  of  these  topis  can 
be  found  at  [9]  and  [10]. 

A.  Range  resolution  by  CW  pulses 

The  NOAA-CODAR  uses  coherent  CW  pulses  for  range 
resolution.  The  length  of  the  pulse  determines  the  spatial  reso¬ 
lution  in  range,  e.g.  a  pulse  with  a  duration  of  8  ps  corresponds 
to  a  1.2  km  wide  circle  around  the  radar.  The  radio  bandwidth 
B  required  for  this  resolution  is  B  =  1/(8  ps)  =  125  kHz. 
However,  due  to  time  multiplexing  ranges  and  antennas  (cf. 
figure  1)  the  effective  sampling  rate  for  a  specific  range  cell 
and  antenna  is  much  less:  488.3  Hz  in  this  case.  Thus,  aliasing 
can  not  be  avoided  by  a  matching  analog  low  pass  filter.  As  a 
consequence  noise  and  interference  signals  are  aliased  back. 
Adding  up  these  noise  bands  reduces  the  signal-to-noise  ratio 
by  24  dB.  This  is  one  of  the  reasons  why  the  NOAA-CODAR 
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like  all  pulsed  radar  systems  uses  high  power  pulses,  e.g.  5  kW 
with  an  average  of  100  W. 


t  =  0  t  =  51  2us 


Fig.  1 .  The  timing  used  within  the  NOAA-CODAR:  64  ranges  and  4  antennas 
are  multiplexed  with  time. 


Advantages  of  the  pulse  technique  are  its  simplicity  and  the 
fact,  that  an  echo  is  always  sampled  into  the  correct  range  cell, 
even  if  it  is  scattered  from  a  target  producing  a  Doppler  shift 
above  Nyquist,  e.g.  a  ship  at  high  speed.  In  this  case  the  echo 
is  aliased  to  the  wrong  Doppler  shift.  Disadvantages  are  the 
high  peak  power  needed  and  the  bad  compatibility  with  other 
radio  services  [13]. 

B.  Range  resolution  by  FMCW 

Performing  the  range  resolution  in  frequency  domain,  solves 
some  of  the  problems  described  above.  By  using  a  conti¬ 
nuously  transmitted  signal  which  is  linearily  increasing  in 
frequency  with  time,  an  echo  at  a  time  delay  At  will  appear  at 
a  constant  frequency  offset  of  A /  (cf.  figure  2).  The  frequency 
of  the  chirp  starts  at  /o  and  increases  to  /0  4  B  during  the 
time  T,  B  being  the  bandwidth  of  the  chirp  and  c  the  speed 
of  light.  This  process  maps  a  target  at  the  distance  r  to  a 


frequency  offset 


A  f  = 


B2r 

t7 


The  received  signal  is  mixed  (multiplied)  with  the  actual 
transmitted  signal  providing  the  frequency  offset.  This  signal 
is  transformed  by  a  Fast  Fourier  Transform  (FFT)  to  resolve 
ranges.  The  spectral  lines  of  this  range  resolving  FFT  are 
quantized  to 


A  q.  _  1  B  2 r quant 

J  quant  rj-t  rj-i  ^ 

In  terms  of  range  cell  depth,  this  frequency  quantization 
corresponds  to 

c 

Tquant  —  %B  * 


A  range  resolution  rquant  =  1.2  km  again  requires  a  band¬ 
width  of  B  =  125  kHz.  To  avoid  range  smearing,  a  windowing 
function  has  to  be  applied  prior  to  each  FFT  [10].  By  tracking 
the  phase  of  consecutive  chirps,  the  Dopper  spectra  within  the 
ranges  can  be  resolved. 

The  advantages  of  FMCW  modulation  are  the  low  conti¬ 
nuously  transmitted  power  (30  W  in  the  case  of  WERA),  and 
a  better  compatibility  with  other  radio  services  as  compared 
to  CW  pulses.  Disadvantages  are  the  required  high  dynamic 
range  of  the  receiver,  which  has  to  handle  the  strong  signal 
from  the  direct  path  (transmit  antenna  to  receive  antenna)  and 
the  weak  signals  from  far  ranges  simultaneously  as  well  as 
a  more  complex  way  of  signal  processing.  Also,  an  echo 
generating  a  Doppler  shift  above  Nyquist  frequency  will  be 
shifted  to  a  wrong  range  cell.  Modern  chips  and  software 
techniques  meanwhile  solved  most  of  these  difficulties. 


C.  Azimuthal  resolution  by  Direction  Finding 

Two  different  techniques  to  resolve  the  incident  angle  of  a 
sea  echo  are  Direction  Finding  and  Beam  Forming.  Direction 
Finding  makes  use  of  the  amplitude  and/or  phase  characteri¬ 
stics  of  a  signal  at  multiple  receive  antennas.  Figure  3  shows 
the  NOAA-CODAR  approach  using  4  antennas  in  a  square 


Fig.  2.  A  linear  frequency  chirp  starting  at  /o  with  the  bandwidth  B  and 
a  duration  T.  After  reaching  the  end  frequency  /o  +  B ,  the  chirp’s  phase 
contiguously  starts  again  at  /q. 


Fig.  3.  A  sea  echo  arriving  from  the  upper  left  which  causes  a  phase 
difference  between  the  antennas. 
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with  A/2  or  slightly  less  diagonal  distance.  An  echo  directed 
from  the  upper  left  in  the  figure  arrives  at  the  4  antennas  with  a 
short  time  delay,  which  can  be  measured  as  a  phase  difference. 
Signals  coming  from  different  directions  are  identified  by  a 
changing  set  of  phase  differences. 

If  phase  information  is  used,  the  incidence  angle  of  a  signal 
can  be  resolved  by  a  least-squares-fit  of  the  measured  phase 
differences  to  the  calibrated  phases: 

4 

e(0)  =  5>*-<M0))2 

i=l 

e(0)  ==>  Minimum 

with  ip*  being  the  measured  set  of  phase  difference,  ipi(O) 
being  the  calibrated  values,  and  0  being  the  incidence  angle. 
The  sum  e(0)  goes  to  a  minimum  at  the  most  propable 
incidence  angle. 

In  absense  of  noise,  the  measured  set  of  phase  differences 
and  the  calibrated  values  are  identical,  giving  the  exact 
solution.  With  increasing  noise,  the  fitted  incidence  angle 
gets  more  and  more  uncertain,  i.e.  the  azimuthal  resolution 
is  a  function  of  signal-to-noise.  If  the  antenna  patterns  are 
distorted,  i.e.  there  are  deviations  from  the  theorectical  func¬ 
tions  of  phase  differences,  additional  systematic  errors  can  be 
introduced,  even  if  the  distorted  patterns  are  known  and  have 
been  taken  into  account. 

An  algorithm  based  on  amplitudes  requires  a  unique  relati¬ 
onship  between  azimuthal  angle  and  amplitude.  For  example, 
two  loop  antennas  installed  perpenticular  to  each  other  like 
with  D.  Barrick’s  SeaSonde  [11]  could  be  used.  Recently,  an 
algorithm  based  on  MUltible  Signal  Characterization  (MU¬ 
SIC)  [12]  has  been  implementedby  D.  Barrick  for  the  Sea- 
Sonde.  All  these  crossed-loop  algorithms  are  patented  by 
“CODAR  Ocean  Sensors”,  but  their  performance  has  not  been 
evaluated  by  the  authors  of  this  paper. 

Non  of  the  published  direction  finding  techniques  is  capable 
of  resolving  the  second-order  sea  echo  spectrum,  which  is 
located  around  the  two  first  order  peaks.  Also,  the  algorithms 
run  into  trouble  when,  at  a  given  Doppler  shift,  sea  echos  from 
multible  directions  superpose.  If  the  number  of  superposed 
directions  is  limited  and  the  signal-to-noise  ratio  is  high, 
MUSIC  can  help  to  solve  this  case. 

I).  Azimuthal  resolution  by  Beam  Forming 

If  access  to  the  complete  backscatter  Doppler  spectrum  from 
a  specific  patch  of  the  sea  surface  is  required,  e.g.  to  get  the 
second  order  returns  for  measuring  ocean  wave  directional 
spectra,  a  linear  array  and  beam  forming  should  be  used.  As 
always,  nothing  is  for  free,  and  the  price  to  pay  in  this  case  is 
the  increased  amount  of  space  needed  to  install  the  linear  array. 
A  16  antenna  array  operated  at  30  MHz  with  A/2  spacing 
between  the  antennas  requires  a  75  m  long  patch  along  the 
shore. 

By  shifting  the  phase  between  antenna  signals  and  adding 
them  up,  a  beam  can  be  steered  to  about  ±60  degrees  to  a 
line  perpendicular  to  the  array.  Figure  4  shows  the  antenna 


Fig.  4.  The  antenna  pattern  of  a  linear  array  with  16  antennas  at  A/2  spacing 
(right  side).  A  windowing  function  has  to  be  applied  to  reduce  side  lobes  (left 
side). 


pattern  of  a  16  antenna  array.  To  reduce  side  lobes  in  the 
antenna  pattern,  a  windowing  function  has  to  be  multiplied  to 
the  amplitudes  of  the  antennas  before  the  signals  are  shifted 
in  phase  and  added  up  to  steer  the  beam. 

A  linear  array  has  a  large  aperture  and  is  quite  insensitive 
against  distortions  due  to  the  environment.  In  many  cases 
the  beam  forming  gives  more  reliable  results  than  direction 
finding,  and  for  access  to  the  full  backscatter  Doppler  spectrum 
it  is  mandatory  in  any  case. 

III.  The  HF  radar  WERA 

The  University  of  Hamburg  WERA  has  been  initially 
developed  in  1996  within  the  EU  funded  Surface  Current 
And  Wave  Variability  Experiment  (SCAWVEX)  project.  One 
main  difference  to  the  NOAA-CODAR  is  the  use  of  FMCW 
modulation  for  range  resolution  instead  of  pulses.  Azimuthal 
resolution  can  be  achieved  by  software  (beam  forming  or 
direction  finding),  depending  on  the  actual  receive  antenna 
configuration  and  software  modules  applied.  WERA  consists 
of  several  hardware  components  which  can  be  put  together 
in  a  modular  way,  i.e.  it  is  possible  to  start  with  a  small  4- 
antenna  system  and  direction  finding  and  upgrade  to  a  16- 
antenna  system  later  on. 


Fig.  5.  The  block  diagram  of  the  WERA  system. 
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Fig.  6.  Radial  components  of  a  surface  current  field  measured  by  WERA  at 
27.65  MHz. 


Fig.  7.  Radial  components  of  a  surface  current  field  measured  by  WERA  at 
16.05  MHz. 


A.  WERA  hardware  components 

The  heart  of  the  system  is  a  low-noise  180  MHz  crystal 
oscillator,  which  is  used  as  a  master  clock  for  all  frequency 
generation  and  sampling  (cf.  figure  5).  This  concept  makes  the 
whole  system  strictly  synchronized  and  coherent.  There  are 
two  independent,  synchronized  DDS  chirp  generators,  which 
produce  the  transmit  and  receive  chirps  at  the  desired  working 
frequency.  In  this  way  constant  frequency  offsets  between 
transmit  and  receive  chirps  can  be  programmed.  The  transmit 
signal  is  amplified  to  an  output  power  of  30  W. 

There  is  one  receiver  module  for  each  receive  antenna.  Each 
receiver  includes  an  HF  band  pass  filter,  an  I/Q  demodulator, 
and  lowpass  filters  to  avoid  aliasing.  High  dynamic  range 
components  have  been  used  to  simultaneously  handle  the 
strong  signal  from  the  direct  path  and  the  weak  sea  echos 
from  far  ranges.  Amplitude  variations  and  phase  shifts  between 
the  antennas  and  receivers  are  measured  and  compensated  by 
software. 

Each  A/D  converter  module  handles  the  demodulated  I/Q 
signals  of  four  receivers.  The  software  automatically  scans 
and  configures  the  correct  number  of  installed  receiver  chan¬ 
nels  in  steps  of  four.  Besides  programming  the  DDSes  and 
controlling  the  A/D  converters,  WERA’s  real-time  computer 
on-line  performs  the  range  resolving  FFTs.  The  real-time 
computer  is  connected  to  a  Linux  workstation  which  controls 
the  measurement  cycles  and  stores  and  processes  the  measured 
data. 

B.  WERA  software  components 

Basically  there  are  three  software  packages  to  run  WERA: 
The  real-time  software  on  the  WERA  itself,  the  WERA  control 
and  calibration  software  which  runs  on  the  Linux  workstation. 


and  the  data  processing  software  which  runs  on  the  Linux 
workstation,  too. 

There  is  a  human  interface  “WeraDesk”  to  control  and  set¬ 
up  WERA  measurement  cycles.  This  software  provides  a  web- 
based  interface  which  can  be  accessed  remotely  through  the 
Internet  depending  on  the  configuration  options  selected  for 
the  Apache  web  server  which  runs  on  the  Linux  workstation. 

Besides  controlling  the  measurement  cycle-  and  repetition 
times,  the  following  radar  characteriscicts  can  be  set:  The 
range  cell  depth  can  be  set  to  3.0  km. .  .0.3  km,  the  number  of 
range  cells  to  32. .  .256,  the  samples  per  data  run  to  64. .  .4096, 
and  the  sample  rate  to  0.173333  s  (good  to  resolve  Bragg 
lines  at  68  MHz). .  .0.520000  s  (Bragg  lines  at  7.5  MHz).  The 
working  frequency  can  be  set  to  anywhere  between  7.5  MHz 
and  68  MHz,  but  the  correct  filters  have  to  be  installed  inside 
the  receiver  modules  and  the  transmit  power  amplifier.  Also 
the  antennas  must  be  designed  for  the  configured  working 
frequency. 

The  data  processing  software  performs  azimuthal  resolution 
(there  are  algorithms  for  direction  finding  and  beam  forming 
as  described  before)  and  calculates  current  maps  or  ocean 
wave  information.  Other  software  modules  can  visualize  the 
measurement  or  calibration  results.  If  there  are  data  links  to 
the  WERA  sites  available,  a  real-time  monitoring  system  can 
be  set  up. 

C.  WERA  measurement  examples 

Up  to  now,  WERA  has  been  operated  at  29.85  MHz, 
27.65  MHz,  and  16.05  MHz.  In  the  very  near  future, 
12.477  MHz  will  be  used  during  an  experiment.  As  expected, 
the  attenuation  of  the  backscattered  signals  increases  with 
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Fig.  8.  The  main  components  of  an  Operational  Forecasting  System  and 
some  examples  of  users  of  such  a  system. 

frequency,  giving  a  working  range  of  up  to  60  km  (cf.  figure  6) 
at  27  MHz  and  up  to  150  km  at  16  MHz  (cf.  figure  7).  At 
12.5  MHz  we  expect  working  ranges  up  to  200  km.  The  range 
resolution  of  WERA  was  possible  to  be  set  to  values  as  low 
as  300  m  at  27  MHz  and  1.2  km  at  16  MHz. 

IV.  An  Ocean  Monitoring  System 

In  the  frame  of  GOOS2  and  its  European  component  Eu- 
roGOOS,  operational  forecasting  of  current  and  wave  fields 
in  coastal  regions  got  more  and  more  important  in  the  last 
decades,  both  for  coastal  management  and  for  security  aspects. 
A  good  overview  on  EuroGOOS  activities  can  be  found  at  [5]. 
One  of  the  key  components  in  this  context  are  high-resolution 
numerical  models,  which  however  require  accurate  forcing  and 
handling  of  the  boundary  conditions.  HF  radar  remote  sensed 
current  and  wave  fields  can  significantly  increase  the  data 
quality  of  the  model  products  through  data  assimilation.  In 
some  cases,  when  the  oceanographic  processes  induce  high 
local  variability,  such  as  mesoscale  eddies  and  fronts,  this 
approach  might  be  the  only  way  to  provide  reliable  now-  and 
forecasts. 

The  general  structure  of  a  monitoring  system  is  given  in 
figure  8.  The  aim  is  to  provide  accurate  on-line  access  to  the 
actual  situation,  where  on-line  means  a  delay  of  one  hour  at 
maximum,  and  to  provide  forecasts.  This  could  be  achieved 
by  a  synergy  of  observations  and  numerical  models,  i.e.  by 
linking  radar  based  measurements  to  fine-resolution  models 
by  data  assimilation.  The  measured  data  are  required  to  force 
the  model  close  to  nature,  whereas  the  model  is  needed  for 
interpolation  and  forecasts. 

2Global  Ocean  Observing  System 


Fig.  9.  EuroROSE  Fedje  experiment:  The  surface  current  field  measured  by 
HF  radar  at  the  3rd  March  2000  5:15  UTC.  A  strong  meandering  front  can 
be  seen  west  off  the  Lyngoy  site. 

On  the  measurement  side,  there  is  e.g.  the  already  establis¬ 
hed  network  of  weather  stations.  This  network  is  delivering 
data  to  the  meteorological  models,  which  provide  the  results  to 
the  weather  services.  The  quatilty  of  these  services  especially 
regarding  the  2  to  3  day  forecasts  has  significantly  increased 
during  the  last  10  years  due  to  the  synergy  of  observations 
and  numerical  models.  One  of  the  aims  of  GOOS  is  to  extend 
the  monitoring  system  from  the  atmosphere  to  the  ocean  by 
including  the  observations  of  buoys,  moorings,  drifters,  and 
remote  sensing  techniques  like  satellites  and  ground-based 
HF  radars  to  the  numerical  models  of  the  ocean  by  data 
assimilation. 

Figure  8  also  lists  some  end  users  of  such  a  monitoring 
service,  e.g.  Vessel  Traffic  Services  (VTS),  port  authori¬ 
ties,  captains  and  pilots  on  sailing  ships,  search  and  rescue 
operations,  and  environment  authorities.  An  example  of  a 
monitoring  system  especially  designed  for  VTS  operations 
has  been  demonstrated  within  the  European  Radar  Ocean 
Sensing  (EuroROSE)[6]  project,  where  the  safe  navigation  of 
oil  tankers  between  the  Norwegian  islands  has  been  supported. 

During  the  EuroROSE  experiments,  two  WERA  HF  radars 
have  been  installed,  north  and  south  of  the  entry.  The  distance 
between  the  WERAs  was  about  13  km.  Figure  9  shows  a 
surface  current  map  measured  by  WERA. 

On  the  model  side,  a  nested  approach  has  been  used.  There 
was  a  three  step  model  chain  for  currents:  The  outer  model 
covers  the  North  Atlantic  and  the  Norwegian  Sea  with  a 
resolution  of  about  20  km.  This  model  delivers  boundary 
conditions  to  an  intermediate  model  (4  km  resolution)  of  the 
coastal  waters  of  southern  Norway,  which  in  turn  provides 
boundary  conditions  to  the  high-resolution  EuroROSE  model 
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Fig.  10.  EuroROSE  Fedje  experiment:  The  surface  current  field  calculated 
by  the  model  after  data  assimilation  of  the  HF  radar  current  field  shown  in 
figure  2. 

of  the  target  area  (60  km  *  60  km). 

The  numerical  models  have  been  operated  by  the  Norwe¬ 
gian  Meteorological  Institute  (met.no),  the  data  assimilation 
technique  has  been  developed  by  the  Norwegian  Nansen 
Environmental  and  Remote  Sensing  Center  (NERSC).  The 
model  and  data  assimilation  technique  is  described  by  Breivik 
and  Saetra  [2], 

Figure  10  shows  a  current  held  delivered  by  the  model/data 
assimilation  system.  The  oceanographic  front  can  still  be 
seen,  although  it  appears  to  be  smoothed.  The  model  results 
represent  the  top  10  m  of  the  sea  surface  (this  is  affecting  the 
navigating  ships),  while  figure  9  shows  the  measured  current 
velocity  at  the  very  top  0.5  m. 

To  get  an  estimate  of  the  performance  of  the  system,  now¬ 
and  forecasts  have  been  compared  to  actual  measurements. 
As  expected,  the  nowcast  and  the  measurement  show  an 
rms  error  as  low  as  10  cm/s  for  a  position  in  the  center 
of  the  measurement  area.  When  comparing  the  2-  to  6-hour 
forecasts  with  the  measurements  taken  at  that  time,  the  rms 
error  increases  to  20  cm/s.  Breivik  and  Saetra,  2001,  present 
scatter  plots  and  correlations  for  the  different  forecast  steps. 
The  correlation  factor  is  0.89  for  the  nowcast,  0.85  for  the 
2-hour  forecast,  0.77  for  the  4-hour  forecast,  0.63  for  the 
6-hour  forecast,  and  0.27  for  a  free  running  model  without 
data  assimilation.  The  last  number  shows  the  importance  of 
measured  data  to  be  included  into  the  monitoring  system. 

V.  Conclusions 

The  University  of  Hamburg  HF  radar  WERA  has  been 
developed  in  1996  and  is  based  on  20  years  of  experience  with 
development  and  application  of  HF  radars.  It  offeres  a  modular 


design  which  can  easily  be  adopted  to  the  requirements  of 
different  applications.  When  operated  at  16  MHz,  a  working 
range  of  up  to  150  km  for  current  measurements  can  be 
obtained.  HF  radars  play  an  increasingly  important  role  within 
monitoring  systems  of  the  ocean  to  provide  accurate  now-  and 
forecasts  to  customers  like  Vessel  Traffic  services. 

Acknowledgment 

This  work  has  been  partly  supported  by  the  European  Com¬ 
mission,  DG  XII,  within  the  Mast-2  programme,  project  CT94- 
0103,  SCAWVEX  (Surface  Current  And  Wave  Variability 
Experiment)  and  the  Mast-3  programme,  project  CT98-0168, 
EuroROSE  (European  Radar  Ocean  Sensing).  The  authors 
are  grateful  to  the  other  partners  of  the  EuroROSE  group: 
H.  Gunther,  W.  Rosenthal,  and  M.  Stawarz  (GKSS,  Germany), 
J.  Guddal  and  M.  Reistad  (DNMI,  Norway),  G.  Evensen  and 
0.  Breivik  (Nansen  Center,  Norway),  L.  R.  Wyatt  and  J.  Green 
(University  of  Sheffield,  UK),  J.  C.  Nieto  Borge  (Puertos 
del  Estado,  Spain),  and  K.  Reichert  (Ocean WaveS,  Germany). 
We  also  wish  to  thank  the  former  members  of  the  HF-radar 
group,  Georg  Antonischki  and  Florian  Schirmer,  as  well  as 
our  technician  M.  Hamann  for  their  support. 

References 

[1]  D.  E.  Barrick,  M.  W.  Evans,  B.  L.  Weber,  “Ocean  surface  current  mapped 
by  radar”,  Science ,  vol.  198,  pp.  138-144,  1977. 

[2]  0.  Breivik  and  0.  Saetra,  “Real  time  assimilation  of  HF  radar  currents 
into  a  coastal  ocean  model”,  Journal  of  Marine  Systems  (JMS),v ol.  28, 
no.  3-4,  pp.  161. ..182,  2001. 

[3]  H.-H.  Essen,  M.  M.  Janopaul,  E.  Mittelstaedt,  J.  Backhaus,  “Surface 
currents  in  the  German  Bight  measured  by  backscattered  radio  waves  -  a 
comparison  with  conventional  measurements  and  model  results”,  North 
Sea  Dynamics  (ed.  by  Siindermann  and  Lenz),  pp.  159-165,  Springer 
Verlag,  1983. 

[4]  H.-H.  Essen,  K.-W.  Gurgel,  F.  Schirmer,  “Tidal  and  wind-driven  parts 
of  surface  currents  as  measured  by  radar”,  Deutsche  Hydro graphische 
Zeitschrift,  36,  pp.  81-96,  1983. 

[5]  N.  C.  Flemming,  Editor-in-Chief,  “Operational  Oceanography  -  Imple¬ 
mentation  at  the  European  and  Regional  Scales”,  Proceedings  of  the 
second  international  Conference  on  EuroGOOS,  ISBN:  0  444  50391  9, 
2002. 

[6]  H.  Gunther,  G.  Evensen,  J.  Guddal,  K.-W.  Gurgel,  J.C.  Nieto  Borge, 
L.R.  Wyatt,  “European  Radar  Ocean  Sensing”,  In  the  proceedings  of 
the  EurOCEAN  2000:  The  European  Conference  on  Marine  Science 
and  Ocean  Technology  of  the  European  Commission,  29.  August  - 
2.  September  2000,  Hamburg,  Germany,  pp.  443-448,  2000. 

[7]  K.-W.  Gurgel,  “Shipborne  measurement  of  surface  current  fields  by  HF 
radar”,  L  'Onde  Electrique,  vol.  74,  no.  5,  pp.  54-59,  1994. 

[8]  K.-W.  Gurgel,  G.  Antonischki,  H.-H.  Essen,  T.  Schlick,  “Wellen  radar 
(WERA),  a  new  ground-wave  based  HF  radar  for  ocean  remote  sensing”, 
Coastal  Engineering ,  vol.  37,  pp.  219-234,  1999. 

[9]  K.-W.  Gurgel,  H.-H.  Essen,  S.  R  Kingsley,  “HF  radars:  Physical  limitati¬ 
ons  and  recent  developments”,  Coastal  Engineering ,  vol.  37,  pp.  201-218, 
1999. 

[10]  K.-W.  Gurgel,  H.-H.  Essen,  and  T.  Schlick,  “The  University  of  Hamburg 
WERA  HF  Radar  -  Theory  and  Solutions”,  First  International  Radiowave 
Oceanography  Workshop,  ROW  2001,  Proceedings,  2001. 

[11]  B.  J.  Lipa  and  D.  E.  Barrick,  “Least-squares  methodes  for  the  extraction 
of  surface  currents  from  CODAR  crossed-loop  data:  Application  at 
ARSLOE”,  IEEE  Journal  of  Oceanic  Engineering,  vol.  OE-8,  pp.  226- 
253, 1983. 

[12]  R.  O.  Schmidt,  “Multible  emitter  location  and  signal  parameter  estima¬ 
tion”,  IEEE  Trans.  Antennas  Propagat.,  vol.  AP-34,  pp.  276-280,  1986. 

[13]  E.  D.  R.  Shearman  and  R.  R.  Unsal,  “Compatibility  of  High  Frequency 
Radar  Remote  Sensing  with  Communications”,  IEE  Conference  Publica¬ 
tion  No.  188,  “Radio  Spectrum  Conservation  Techniques”,  pp.  103-107, 
IEE,  London,  July  1980. 


705 


